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We enlarge the classes of infiaton and quintessence fields by generalizing the pseudo 
Nambu-Goldstone boson potential by means of elliptic Jacobian functions, which are 
characterized by a parameter k. We use such a generalization to implement an inflation¬ 
ary era and a late acceleration of the universe. As an inflationary model the Jacobian 
generalization leads us to a number of e-foldings and a primordial spectrum of perturba¬ 
tions compatible with the Planck collaboration 2015. As a quintessence model, a study 
of the evolution of its equation of state (EoS) and its w'-w plane helps us to classify it 
as a thawing model. This allows us to consider analytical approximations for the EoS 
recently discovered for thawing quintessence. By using JLA Supernovae la and Hubble 
parameter H[z) data sets, we perform an observational analysis of the viability of the 
model as quintessence. 

Keywords: Dark energy; observational cosmology. 
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1. Introduction 

Scalar fields are widely used to endeavor to explain some aspects of the modern 
cosmology, in particular, the inflationary epoch in early times and the recent accel¬ 
eration of the universe. Mechanisms for inflation and recent acceleration are not so 
different. The main difference is that non-relativistic matter cannot be ignored in 
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recent times. Moreover, the energy scale of the quintessence potentials needs to be 
much smaller than the inflationary ones. 

Inflation is a period preceding the standard cosmology where the evolution of 
the universe can be described by a scalar field, the inflaton. It was proposed to solve 
the horizon and the flatness problems of the standard model, as well as provide the 
mechanism to generate the seeds for the formation of structures and anisotropies of 
the cosmic microwave background (CMB). During the inflationary period, the uni¬ 
verse accelerates its expansion to an almost exponential rate, cooling and smoothing 
the spatial structure. Furthermore, inflation predicts fluctuations with an almost 
scale invariant spectrum which has been confirmed by recent CMB data.®^ To con¬ 
nect with the standard cosmological periods dominated by radiation and matter, it 
is necessary to add a period of reheating. The usual mechanism can also be preceded 
by a more efficient phase known as preheating.^^lt^ 

On the other hand, several cosmological experiments have demostrated the ex¬ 
istence of a phase of accelerated expansion in current times.^^^^^^^ This leads us to 
one of the most challenging problems in cosmology: the nature of dark energy. So 
far, the most sucessful dark energy candidate is the cosmological constant A. How¬ 
ever, this model apparently has shortcomings such as discrepancy between the value 
of the vacuum energy obtained through observations and its estimated theoretical 
value, or even the coincidence problem.l^^ti^ Thus, in recent years, alternatives to 
solve these problems have been sought. Alternatives to the ACDM model can be 
divided into two groups. In the first group, general relativity is modified to obtain 
late acceleration without exotic components (see for example Refs. fT^lfTSl) . The sec¬ 
ond group is related to the existence of an exotic fluid with negative pressure such 
as Chaplygin gas, /c-essence, quintessence and othersin all these alternatives 
the dark energy equation of state (EoS) changes dynamically with time differently 
from the ACDM, the EoS of which is constant and equals —I. In the context of 
quintessence, a canonical scalar field (j) minimally coupled to gravity is used as a 
dark energy model. The field 4> varying slowly along a potential E((/>) can lead to 
observational results very similar to the cosmological constant (for a review see for 
example Refs. [Ml[2^ . 

One of the most interesting models for inflation and dark energy is an ultra-light 
pseudo Nambu-Coldstone boson (pNCB) which is still relaxing to its vacuum state. 
From the viewpoint of quantum field theory, pNCB models are the simplest way 
to have a naturally ultra-low mass, spin-0 particles and hence, perhaps the most 
natural candidate for a minimally coupled scalar field presently existing. The pNCB 
effective potential is given by 



( 1 ) 


where / ^ nip ~ 10^® CeV is the Planck scale representing the spontaneous breaking 
scale, and /i is related to the explicit global symmetry breaking scale. Note that 
/ —>■ 00 corresponds to the exact shift symmetric constant potential. pNCBs were 
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first proposed in the context of natural inflationf^ and then they were subsequently 
extended for the case of dark energy!^ Also, a semiclassical calculation of particle 
production by a Nambu-Goldstone boson is performed in Ref. |28l 

On the other hand, analytical integration of the motion equations, in context of 
both inflation and quintessence, admits as solutions potentials related to hyperbolic 
functions.l^SHlll In fact, some of these potentials can be considered as hyperbolic ex¬ 
tensions of O- Additionally, elliptic functions have been considered to study some 
inflationary aspect^^ and it was recently shown, in context of Hamilton-Jacobi 
method, that solutions related to the inflationary potential can be expressed in 
terms of elliptic functions.^^^ Furthermore, in recent times, Friedman-Robertson- 
Walker universe (FRW) containing a mixture of perfect fluids and a quintessence 
escalar field, admits as solutions for the potential, some especial cases of elliptic 
functions.^ Within this context and keeping in mind that trigonometric and hy¬ 
perbolic functions are particular cases of Jacobi’s elliptic functions, we propose, in 
some sense, a unified treatment of all of these potentials. For this reason we gen¬ 
eralize Eq. ([T]) in an elliptic version by means of Jacobi function^, and we use it 
to implement an inflationary era and a late acceleration of the universe. Such a 
generalization owns a shift symmetry, like the pNGB potential, important to solve 
the two reasons that render challenging to find a natural candidate for quintessence 
field in particle physics. The first reason is related to the violation of the flatness 
condition (see for example Ref. l45ll . while the second is related to the fact that an 
ultra-light field would carry a fifth force, typically of gravitational strength (see for 
example Refs. I46I47|I . In this work we shall limit ourselves to the case of the FRW 
homogeneous and isotropic metric with null curvature. 

This paper is organized as follows: in Section[^we present the generalized poten¬ 
tial which can be expressed in terms of elliptic functions. In Section[3]an inflationary 
model based on the generalized potential is presented. We provide the general so¬ 
lutions and the expressions for the relevant cosmological parameters and then we 
study the viability of the model in the context of recent data from the Planck col¬ 
laboration 2015.^ In Section [4] we study the generalized potential as a dark energy 
model, i.e., as a quintessence model. We perform the study of the evolution of its 
EoS and of the phase space w' — w . Finally, we perform observational tests to in¬ 
vestigate the viability of quintessence models based on the Jacobian generalization. 
In section [5] we present our conclusions and final remarks. 


^Jacobi’s elliptical functions are not strange in physics and engineering. They appear in a variety 
of problems. For example, they appear when the equation of motion of the bead on a rotating hoop 


is solved,!^ or in the study of the nonlinear Schrodinger equation in the context of Bose-Einstein 
condensates,or when population growth dynamics is studied,!^ as well as in other fields of 

physics 
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2. Jacobian pseudo Nambu—Goldstone boson 

Recalling that trigonometric functions are special cases of elliptic functions, we 
proceed to generalize the effective potential © in a simple way: 


R(<^, k) = / 


1 + cn 



( 2 ) 


where cn(a:) = cn(a::| k) is the Jacobi’s elliptic cosine function, and k is the modulus of 
the elliptic function. Note that the original pNGB potential is recovered by making 
k = 0. We will refer to potential ([2]) as Jacobian pseudo Nambu-Goldstone boson 
(JpNGB) effective potential. In FIG. [I] the effective potential (I2|) is plotted for 
different values of the modulus k and in the Appendix some properties of the Jacobi’s 
elliptic function are presented. Without loss of generality, we will consider that /i 
and / possess the same physical meaning as in the pNGB scheme. 



Fig. 1. Plot of the Jacobian pseudo Nambu-Goldstone boson effective potential V as a function of 
the dimensionless scalar field cp = (f>l f for different values of the modulus k. Clearly, the hyperbolic 
limit k = 1 does not allow inflation to end because the infiaton rolls eternally. 


Alternatively, we mention that Higaki & TakahashP^ have recently studied the 
case with a minus sign in the potential ([2]) calling this model the Cnoidal inflation. 
Specihcally, they have shown that this model predicts values of the spectral index 
and the tensor-to-scalar ratio that interpolates from natural inflation to exponential 
inflation such as R^-or Higgs inflation and brane inflation. In the next section we 
will use potential ([2]) to implement an inflationary epoch in the early universe and 
we will study its concordance with observational data of the Planck satellite. 
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3. JpNGB as a inflaton field 

In the inflationary regime, the dynamical evolution of the JpNGB field is described 
by the Klein-Gordon equation 

^+i3H + T)^ + V^ = 0, (3) 

where F is the decay width of the inflation, the dot indicates the derivative with 
respect to time, H = a/a and represents the derivative of the potential with 
respect to (j>. Note that in the temperature range IJ- T < f, the JpNGB potential 
is dynamically irrelevant because the term is negligible when compared to the 
Hubble damping term. Therefore, in this temperature range, aside from the smooth¬ 
ing of spatial gradients in (j), the field does not evolve. For temperatures T < f, the 
global symmetry is spontaneously broken, and the field (j) describes the phase de¬ 
gree of freedom around the bottom of the potential. Since (j) thermally decouples 
at a temperature T ^ f we assume it is initially laid down at random 

between 0 and Pf in different causally connected regions, where P is the period 
of the elliptic function. Finally, at temperatures T < /i, in regions of the universe 
with (j) initially near the top of the potential, the field starts to roll slowly down 
the hill toward the minimum. In those regions, the energy density of the universe is 
quickly dominated by the vacuum contribution [V{4>) ~ > Prad ^ P], and the 

universe expands exponentially. The slow-roll (SR) regime occurs when the motion 
of the field is overdamped, so ^ <C 3iL(^ and F ^ 110, and therefore two conditions 
are met: 




( 4 ) 


and 


vm 


v{4>) 


< 


\/487r 


( 5 ) 


which, using Eq. ([5]) and defining the dimensionless scalar field p = cj)/f, leads to 


'2cn(^ |dn^ — fc^sn^ ^ \/487r/ 


and 


\ + cap 


|dn(^sn(^| ^ \/487r/ 




( 6 ) 


( 7 ) 


l + CTip nip 

respectively (see definitions and relations of Jacobi’s elliptic functions in the Ap¬ 
pendix). Obviously, the SR regime ends when one of inequalities ([6|) or ([7|) is vio¬ 
lated. This fact occurs when (j) reaches a value (pe which is evaluated for / = rup and 
f = mp/V^ in FIG.H Note the consistency between our results and those found 
previously by Freese, Frieman & OlintcP^ using k=0, i. e. ip\ 




= 1.9 for / = mpl\/2^'K. 


= 2.98 for / = rUp and 


^For warm inflation, the inflaton motion is also overdamped but under condition F > /t. 
















July 8, 2016 0:23 WSPC/INSTRUCTION FILE JpNGB2 


6 W. S. Hipolito-Ricaldi and J. R. Villanueva 



Fig. 2. Plots for the condition m in the (p—k plane for two values of /. The validity of inequality 
0 is represented by the grey region of each graph and the boundary corresponds to the value 
of dimensionless scalar field p at the end of infiation, pe- Also, the value reported by Freese et 
alfm for fc = 0, is represented by a dot. Top panel: / = mp and p^ = 2.98; Bottom panel: 
f = mpl^2A.'K and p^ = 1.9. 
























July 8, 2016 


0:23 


WSPC/INSTRUCTION FILE JpNGB2 


A Jacobian generalization of the pseudo Nambu-Goldstone boson potential 7 



Fig. 3. Numerical solution for the initial value ipi and final value of the inflaton as a function 
of the modulus k using N = 7t) and / = mpl\/2i-K. 


Now let US calculate the values of the inflaton related to the number of e-folds 
required to generate the seeds for the creation of the large scale structures. As is well 
known, the number of e-folds of physical expansion that occur in the inflationary 
stage is given by 


where 




Stt V 

"ip Jf y,<i> 

ds y? / 1 - cn 1 + cn tp ^ 
TOp \dsipe\ 1 -b cn 93e 1 — cn 



( 8 ) 

(9) 


and fc' = 1 — fc is the complementary modulus. Therefore, knowing the value of (jte 
from Eq. o the initial value of the inflaton (j)i can be calculated, imposing a value 
on N, for example, N{(j)i) = 70. In FIG.[3]we show the initial (pi and hnal (pe values 
of inflaton calculated using Eq. dH) and / = mpl'j2ATT as a function of the modulus 
k. Our results indicate that for any value of k, the inflaton with a JpNGB effective 
potential is compatible with 70 e-folds. We have checked for values of between 

60 and 70 and have found that the compatibility is maintained. 

In the inflationary scenario, the quantum fluctuations are relevant because they 
generate two important types of perturbations: density perturbations (arising from 
quantum fluctuations in the scalar held, together with the corresponding scalar 
metric perturbation), and relic gravitational waves (which are tensor metric huctu- 
ations). The former is sensitive to gravitational instability and leads to structures 
formation, while the latter predicts a stochastic background of relic gravitational 
waves. Scalar perturbations produce a power spectrum characterized by the scalar 












July 8, 2016 0:23 WSPC/INSTRUCTION FILE JpNGB2 


8 W. S. Hipolito-Ricaldi and J. R. Villanueva 



0.94 0.96 0.98 1.00 


Primordial tilt rig 


Fig. 4. Plot of the tensor—to-scalar parameter r as a function of the scalar spectral index Us 
for the Jacobian pseudo Nambu-Goldstone Bosons potential for / = mpl%/24w. The marginalized 
joint 68% and 95% confidence level regions using Planck TT + low P, Planck TT, TE, EE + low 
P and Planck 2013 data release with k = 0.05/r Mpc~^ are shown.l^ 


spectral index Ug while tensorial perturbations produce a power spectrum charac¬ 
terized by the gravitational wave spectral index ut- In inflationary models it is 
common to define the tensor-to-scalar amplitude r = where Pt and Pr are 
amplitudes of the tensor and scalar power spectrum respectively. 

The spectral index rig and the tensor-to-scalar amplitude ratio r are written in 
terms of the slow-roll parameters, 


ey = 


ilv = 


ml 


l" ^l 

1 2 9 

an (f sn (f 

IGtt 



(1 + cnpY 

ml / 


] 

cnp{2k sn^ p — Y 

Stt \ 

. V . 

I IOtt/^ 

1 J- cn 


the spectral index and the tensor-to-scalar ratio are given by 

ng ~ 1 -I- 2 Tjy — 6 ey 

SiTip dn^ ip sn^ ip iTip cnp (2 k sn^ p 
Stt/^ (1 +cnp)^ Jtt/^ 
rrip dn^ p sn^ p 


1 -I- cn 


"tri on i/ysn p 
r ~ 16ey = —^ 

(1 + c-o-pY 


( 10 ) 

( 11 ) 


( 12 ) 

(13) 


Therefore, using rig and r, it is possible to obtain a direct comparison of our results 
with the Planck collaboration data sets, which is shown in Figure 01 Based on such 
comparisons, we may say that a description of the inflationary universe model in 
terms of a JpNGB effective potential could quite well accommodate the data recently 
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released by the Planck missionP However, potentials with k ^ 1 are not viable 
because the inflaton rolls eternally. 


4. JpNGB as dark energy 


It is well known that a canonical scalar field called quintessence, varying slowly 
along one potential V{4>) can lead to observational results very similar to the cos¬ 
mological constant and would explain the late acceleration of the universe (for a 
review see for example Ref. [35]). However, two reasons render it challenging to find a 
natural candidate for (j) in particle physics. The first is related to the violation of the 
flatness condition,l2Sl while the second is related to the fact that an ultra-light field 
would carry a fifth force, typically of gravitational strength.SsEll These problems 
can be solved by means of a potential which has a shift symmetry, as the pNGB 
potential.!^ 

The fact that potential ([Sj), which is a generalization of pNGB, is a periodic 
function and therefore has a shift symmetry motivates us regarding the possibility 
of using it in order to construct a viable quintessence model. Thus we exploit its 
dynamic behavior when it is used to explain accelerated expansion in late times. 

We start by considering the general motion equations in the case where the 
universe is dominated by a pressureless matter and a quintessence field. In the 
context of a homogeneous and isotropic Friedmann-Robertson-Walker universe, with 
null curvature and in the case where non-relativistic matter is not interacting with 
the quintessence field we have 


^ + 3H^ + V^ = 0 ( 14 ) 

Pm + iHpm = 0 ( 15 ) 

Pfj, + ZH{p^+p^) = Q (16) 

iH^rn^p = pm -f P(p ■ (17) 


where p^ = + V{(j)) and p,p = — V{(j)) are the density and pressure of 

the quintessence field respectively, and w = p^i/p^, is the EoS of dark energy and 
subindex m indicates matter. From ([T4l) - (II7p equations for w and = p^/pcrit, 
where pcrit = can be written (see for example Refs. jH] jSnj) 


w' = (I 


w) 


—3(1 re) Ay^3(I 

% = -Swci^i - . 


(18) 

(19) 


Additionally, one equation for the auxiliary variable A = —mpV^jV is found 


A' = -^3{i + w)n^x^{r - 1), 


( 20 ) 


where F = V/E| and ' indicates derivative with respect to In a. Thus, given 
a potential V, F is known and the system (IT8l) - (l2n]) is completely determined. For 
example, for the exponential potential F = I and A is a constant and for pNGB 
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Fig. 5. Evolution of EoS it; as a function of scale factor a for different values of the parameter k 
but same initial conditions. The case fc = 0 represents pNGB potential. 

potential F = (A^ — 1)/2A^. For our case F has a more general dependence of A and 
we will find it numerically. 

Recently, it was shown that quintessence models can be approximately sepa¬ 
rated into two classes, depending on the acceleration or deceleration of the field as 
it evolves down its potential towards to a zero minimum.®! Such groups are called 
thawing and freezing quintessence and they behave quite differently. Models where 
the field is nearly frozen away from the minimum of the potential due to large 
Hubble dumping during the early cosmological era and starts to evolve towards the 
minimum recently are called thawing. Models where the field initially rolls towards 
its potential minimum and has only recently gradually slow down are called freez¬ 
ing.® Examples of the first group are the exponential and pNBG potential^^lllMlISj 
and examples of the second group are the inverse power law potentials.®® In the 
next subsection we will show that our generalization belongs to the thawing group. 


4.1. JpNGB as a thawing quintessence 

As noted above, in thawing quintessence models, the field is initially frozen in early 
times and then as H increases the field rolls down the potential. This means that at 
early times the EoS is w ~ —I but grows less negative with time where w' > 0. Since 
we are generalizing the quintessence potential by means of (ED, we are interested in 
knowing if such thawing behavior appears in this particular case. 

In order to search the thawing characteristic we solve numerically the system 
(fT^ - (j20|) with same initial conditions and for different values of fc-parameter. Results 
are shown in Figures © and ([SD. In FIG. © the evolution of the EoS as a function 
of scale factor a is presented. For all cases we have used the same conditions: f /nip = 
0.9 and 4>i/ f = 1.5. As we can see, the same initial conditions lead us to different 
evolutions of EoS for each k. This opens the possibility of discriminating between 
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different members of the family potentials ([2]) using observational data. In all cases 
w is close to —1 for early times but begins to evolve for late times. Of course 
k = 0 represents pNGB potential. We can note also that in the past, models with 
potential @ are indistinguishable from the pNGB model and differences between 
them appear only recently. Additionally, it can be seen that the evolution of w is 
slower for higher values of k. 

In FIG. (IH) we plot the behavior of the model on the dynamical phase space 
w' — w for different values of k and for the same initial conditions. The upper 
shadowed area represents the thawing region I + i(;<w'<3(1 + 1 /;) and the lower 
shadowed area represents the freezing region 0.2//;(I + w) < 3//;(I + w)^ 

The values used in FIG.® are f/mp = 1.1 and 4>i/f = 1.5. Note that in all 
cases, curves are inside the thawing region and we can thus infer that JpNGB has 
thawing quintessence features. We have checked for a wide range of initial conditions 
(0.6 < (j)ilf < 2.5) and fc-parameters (0 < fc < I) and such thawing features are 
maintained. Thus we would expect that thawing features are always maintained. 



Fig. 6. Plane representing the dynamical space phase w' — w and curves for JpNGB effective 
potential for different values of the modulus k. The upper shadowed area is the thawing region 
^ + < 3(1 + w) and the lower shadowed area represents freezing region 0.2tii(l + tii) < u?' < 

3w{l + w). 


Interestingly, for models having freezing and thawing properties an approxima- 
tion of w can be known analytically. In Refs. [48l [58] the equations of motion (ITdll 
were solved in the limit when w is close to —I and an approximated EoS as a func- 
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tion of the scale factor was found when |?/; +1| <C 1. The former condition is ensured 
if conditions of slow roll thawing quintessence 'C 1 and -C ij^ in the 

case where both matter and a scalar field contribnte to the density are satisfied. 
Such approximation gives us (see Refs. l48l[58l for details) 


where 
and it 


u;(a) = -! + (! + x 

( (K-F{a)){F{a) + l)^ + {K + Fia))iF{a)-l)^ T 

\{K- + 1)^ + iK + 

wq is the EoS in the current time and K measures the potential curvature 
is 


K = 


VM^i) 

3 V{^i) 


( 22 ) 


with (f>i being the initial value for (j) and function F(a) has the form 

F(a) = ^l-(r!^J-l)a-3. (23) 

Such approximation was found by means of a Taylor expansion around and 
comparisons with other approximations like Chevallier-Polarski-Lind^^ and Crit- 
tenderPS were performed in Ref. [58l As shown in Refs. 081158] numerical solutions 
are in excellent agreement with the approximate w for several quintessence models 
including pNGB. We made the same comparison between exact and approximated 
EoS for JpNGB quintessence models. Results are presented in FIGl7]for different 
values of k. We have used the following values (i) k = 0.4 f/uip = 2, pi/f = 3 
{K = 1.6); (ii) k = 0.8, f/rup = 2 and pi/f = 4 {K = 1.6) and (iii) k = 0.99, 
f/rup = 2 and pi/f = 6.6 {K = 1.1). In all figures, continuous curves correspond 
to an approximate solution lEa. (l21|) l and dashed curves represent exact numerical 
solutions of the system (IT^ - dSnil . 

In order to evaluate the relative difference between exact and approximate solu¬ 
tions, the relative error (|(5r(;/'u;|) was constructed for the curves in FIG. ([7]). Results 
are presented in FIG. ([8]). Note that in all cases the error is < 0.9%, which is in 
agreement with relative errors for other quintessence models.^ We have examined 
several cases by varying initial conditions, k parameter, etc. and in all of them the 
highest relative error was ~ 1%. As a conclusion, Eq. (l2T]l gives a good approximate 
solution of EoS for our generalized model. Parametrization (EID lead us to wq 
and K as free parameters for EoS. In the next section we will use (l2T]l to study 
observational viability of JpNGB quintessence models and constrain free parame¬ 
ters by using recent observational data. Additionally we will put constraints on the 
fc-parameter. 
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a 


Fig. 7. Evolution of the EoS as a function of a using exact and approximate solutions. Dashed 
curves represent exact numerical solutions while continuous curves represent approximate solu¬ 
tions. 



Fig. 8. Curves with relative error as a function of scale parameter a for curves (i), (ii) and (iii) 
of FIG.llTt. Note that relative error is at most < 0.9%. 

4.2. Observational viability 

Analysis considering Eq. (I^Tl) for thawing models have already been carried out in 
Refs. SHI EU [59l For example in Ref. |6T]the likelihood analysis with the SNIa and 
BAO data was performed. In Ref. [SHI these analyses were updated by using the 
Union 2.1 data setfs3 the CMB shift parameter measured by WMPAtI^ and the 
BAO data of the BOSS experiment.!^ 

In this subsection we present an updated observational analysis by using the 
latest SNIa data (JLA sampl^^ and H{z) data.!^ Also, a first attempt to apply 
Eq. (1^ to study JpNGB models is developed. Our tests are based on y^-statistics, 
which will allow us to explore the space of parameters and constrain the free pa¬ 
rameters 9. 
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Function is given by 

x\0) = Ay^{e)C-^Ay{9), (24) 

where Ay{9) = yi — y(xi\ 9), C is the covariance matrix of data points yi and y{xi\9) 
represents theoretical predictions for a given set of parameters. The best fit is found 
by minimizing the y^-function and the minimum of gives us an indication of the 
quality of the fit. 

First, we considered tests involving the distance modulus of type la Supernovae, 
which is defined by 

yL{z, 0) = 5 logio{dL{z, 9)) + 42.38 — 5 logioh , (25) 

where the luminosity distance is given by 

dt(*,«) = (i+ ^ <*> 

For our observational treatment with SNIa we used the JLA sampl^^ which has 
740 data points. This updates a previous analysis performed iiP^ based on the 
Union2.1 data set.^ Observational data points of the luminosity-distance modulus 
were calculated using relatioiP^ 

Mohs =m*B- {Mb - aXi + pC) , (27) 

where m*g corresponds to the observed peak magnitude in the rest frame B band 
and a, /3 and Mb are nuisance parameters, Xi is related to the time stretching of the 
light curves, and C corrects the color at maximum brightness. In order to calculate 
completely Mohs and its covariance matrix we have followed the steps suggested in 

Ref.lMl 

After that, we performed an analysis by using the compilation of observational 
points for H(z) parameter, which were derived using the differential evolution of pas¬ 
sively evolving galaxies as cosmic chronometers.l^l^ We used the recently updated 
H{z) datcl^and finally we performed an joint analysis using Xtotai ~ XsNiad~XH(z)- 
Considering Eq (1^ in an observational context of fJ,{z,9) and H{z) leads us 
having 4 free parameters (h, wq and K). However, we have marginalized over 
parameter h with the prior 50km s~^Mpc~^ < Hq < 80km s~^Mpc~^ in all of our 
analyses. Thus, we leave all other three parameters free. As discussed irP^ and for 
reasons related to approximation limits of Eq. m, we used a prior 0 < AT < 10. 
Using SNIa alone, our results give us the constraint of = 0.73to;o7: '^o = 
— 1.13^0.17 and K = 0.991® gg (2cr CL). Using H{z) alone our results give us the 
constraints for = 0.72ljj;Q5, wg = —0.80^g;35 and K = O.llgl (2cr CL), while 
the joint analysis gives us H^g = 0.73lQ g 5 , wq = —I.OSlgl® and K = 0.99lQ;gg 
(2tT CL) with a xl ^ l-l- Results of the joint analysis are presented in FIC. ([9|). 
FIC. El (a) shows us the PDF function of the parameter D^g after marginalization 
over Wq and K. FIC. [9] (b) shows us the PDF function for parameter wq after 
marginalization on D^g and K. FIC.|9](c) shows us the PDF function for parameter 
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(a) (b) 




(c) (d) 

Fig. 9. (a) PDF function for parameter after marginalization on wo and K, (b) PDF function 

for parameter wq after marginalization on ^iid K, (c) PDF function of parameter K after 
marginalization on ^ud wq and (d) wq - K plane with contour lines at Icr, 2a and 3cr CL. 


K after marginalization on and wq. Finally FIG. |9] (d) visualizes the wq - K 
plane with contour lines at Icr, 2cr and 3cr confidence levels. Our results here are in 
agreement with results in.l^ It can be seen that there is a very high degeneration 
for parameter K. 

So far, we have considered the parameter K. However, in general, this parameter 
is not adequate to discriminate between different thawing models through the data. 
Since we are interested in knowing if JpNGB models are preferred or ruled out 
by the observations, it is necessary to consider K as & function of parameter k. 
Actually, K is dependent on of also. However, since ipi/f and k are poorly 
constrained by data considered here, we will fix the value of 4>i/f and vary the 
other three parameters. This will allow us to have information about the remaining 
parameters, particularly of the k parameter. 

Thus, considering Eqs. dill), ^ and m we have performed a new analysis 
where free parameters are now ^l^o, wq and k where we have fixed 4>i/f = 1.9. Using 
SNIa alone our results give us constraints for = 0.76lg;Q5, wq = —I.I 6 ^o ;^4 
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and k = 0.62to;62 ( 2 ^ CL). By using H{z) alone our results give us constraints 
for = 0.72fg;Q5, wq = —O.SS^gg® and k = O.Stg;^ (2tT CL), while the joint 
analyses gives us wq = — k = 0.25to;25 CL) with 

xl ^ 1-14. Results of the joint analyses are presented in FIC. dTUl) . FIC.[Tn](a) shows 
us the PDF function for parameter wq after marginalization over and k FIC. 
1101 (b) shows us the PDF function for parameter k after marginalization over 
and Wo • FIC. [TU] (c) shows us the wq - k plane with contour lines at Icr, 2cr and 3 ct 
confidence levels. FICHU] (d) visualizes the - wo plane with contour lines at Itr, 
2ct and 3tT confidence levels. Note that there is a high degeneration in fc-parameter. 

In order to see how our results are affected by <j)i/f, we have varied such a 
parameter in the range 0.6 < (j)i/ f < 2.5 and we redone the analysis. Results indicate 
that best-fit values for and wq vary in ^ 5% as maximum. However, best-fit 
value for k varies ^ 40%. Such high variation is related to the strong degeneration 
in k. 



Fig. 10. (a) PDF function for parameter wq after marginalization over (t>) PDF 

function for parameter k after marginalization over caQ. (c) wq -k plane and (d) ■ ”^0 

plane. In (c) and (d) contour lines at Icr, 2c7 and 3cr CL are presented. 
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5. Conclusions 

In this work an elliptic generalization of the classical pseudo Nambu-Goldstone po¬ 
tential is studied. We have denominated such a generalization the Jacobian pseudo 
Nambu-Goldstone (JpNGB) potential. The JpNGB is used as an inflationary po¬ 
tential and we discuss the general conditions of the parameters to describe an infla¬ 
tionary stage. Additionally, we have calculated the values to the inflaton necessary 
to obtain the number of e-folds of physical expansion required to generate the seeds 
for the creation of the large structures in the ini verse. Also, the relevant quantities 
for scalar and tensor perturbations, i.e., the scalar spectral index rig and the scalar- 
to-tensor ratio r, have been calculated explicitly and then contrasted with recent 
observational data from the Planck collaboration, where marginalized joint 68% and 
95% confident level regions for Planck plus WMAP data for the model AGDM plus 
r were used. These analyses predict that our generalization is able to describe the 
inflationary stage with exception of potentials with /c ^ I. 

Afterward, we have used the JpNGB potential and the construction of a fam¬ 
ily of quintessence models parametrized by parameter k was implemented, where 
pNGB quintessence is the particular case k = 0. We solved equations for generalized 
potential in order to study the evolution of the dynamical EoS and its evolution 
with the scale factor. It was found that the same initial conditions lead us to a dif¬ 
ferent evolution of EoS for each different k, opening the possibility of discriminating 
between them by using observational data. 

As is well known, pGNB is a typical thawing quintessence model. In order to 
see if similar behavior appears for the family of generalized models, we performed a 
study of dynamical phase space w' — w. It was confirmed that such a generalization 
maintains the thawing features of the pNGB for all the possible values of k. This 
behavior allows us to use an analytical approximation of the dynamical w{a) recently 
discovered for thawing models. 

In order to ascertain whether the generalization is compatible with the observa¬ 
tions, we performed an observational analyses using SNIa and H(z) data and the 
analytical solution for w{a) (Eq. (l?T]) b After marginalization over h, the space of 
parameters in this case is K). The best-fit values using SNIa and H{z) 

joint analyses are = O-^SIq q^, wq = — I-051q and K = 0.99^0^9 (2cr GL) 
with sl Xu ^ l-l- These results are in agreement with previous analyses. However, in 
general, the AT-parameter is not adequate to discriminate between different thawing 
models through the data. Since we have been interested in knowing if some general¬ 
ized model is preferred or ruled out by the observations, it is necessary to consider 
AT as a function of k. Thus, after marginalization over h and using a joint analyses 
with SNIa and H{z), we have as best-fit values = O.TStg^gij, wo = —I-06j|jQ;i[J 
and k = 0.25^0.25 CL) with a yj ^ 1.14. These results show us that at 2cr 
GL, JpNGB models with 0 < fc < I are not ruled out by SNIa and H{z) data,i.e., 
quintessence fields with potential (l2|) with the ~ 70% of the content of the dark 
energy and with u> ~ — I are compatible with SNIa and H{z) data. Nevertheless it 
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seems to exist a strong degeneration in the fc-parameter which needs to be broken 
using another data sets and with a detailed study at perturbative level. 

Finally, we can conclude that the generalization considered here is a viable model 
and it works well overall for the inflationary era and for the acceleration era in late 
times. Of course a better and a more detailed study is needed that considers the 
limits of approximation for (1211) and exact numerical solutions as well as other data 
sets like CMB or matter power spectrum. Such points are still under investigation. 
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Appendix A. A brief review of Jacobian elliptic functions 

The Jacobi elliptic functions are set of elliptic functions whose three basic functions 
are called the Jacobi elliptic cosine cn (u, fc), Jacobi elliptic delta dn (u, fc), and Jacobi 
elliptic sine sn (it, fc), where k is known as the elliptic modulus. They arise from the 
inversion of the elliptic integral of the first kincp2HI3 

u{y, k) = u= f 
Jo 


/O V(l-t2)(l-fct2) 

de 


= F{ip,k), 


(A.l) 


Jo \/l — k sin^ 9 

which was studied and solved by Abel and Jacobi. The trigonometric and hyperbolic 
limit are 


sn(ii, 0) = sinit, sn(ii, 1) = tanhit, 
cn(it, 0) = cosit, cn(it, 1) = sechu, 
dn(M,0) = l, dn(it, 1) = sechu. 


(A.2) 

(A.3) 

(A.4) 


The quotients and reciprocal of {snu, cnu, dnu} are designated in Glaisher’s no¬ 
tation by 


nsu = 

1 

csu = 

cnu 

dsu = 

dnu 

(A.5) 

snu’ 

snu ’ 

snu ’ 

nc u = 

1 

scu = 

snu 

dcu = 

dnu 

(A.6) 

cnu’ 

cnu’ 

cnu ’ 

ndu = 

1 

sdu = 

snu 

cd u = 

cnu 

(A.7) 

dnu’ 

dnu’ 

dnu 
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Therefore, in all, we have twelve Jacobian elliptic functions. Finally, some useful 
fundamental relations between Jacobian elliptic functions are 


2,2 1 
sn M J- cn u = 1, 

(A.8) 

dn^u + k sn^u = 1, 

(A.9) 

dn^u — k CT^u = k', 

(A.IO) 

2 1 7 / 2 J 2 

cn u + k sii u = an u. 

(A.ll) 
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